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1.  INTRODUCTION 

In  most  practical  data  gathering  tasks ,  the  removal  of 
unwanted  frequency  components  from  the  data  by  filtering  is  an  essential 
step.  As  the  cost  of  microprocessors  decreases  and  as  their  speed  and 
capabilities  increase,  digital  filters  implemented  by  microprocessors 
become  an  attractive  alternative  to  analogue  filters,  especially  in 
low  frequency  data  gathering  applications. 

Although  there  is  a  great  deal  of  information  available  on 
the  design  of  digital  filter,  it  tends  to  be  in  a  form  which  is  either 
too  limited  in  scope  fl,2l  or  too  tedious  or  difficult  to  read  for 
the  average  reader  ^33,5'~j  . 

This  memorandum  is  an  attempt  to  overcome  this  problem  by 
choosing  a  coononly  used  filter  configuration,  the  digital  Butterworth 
low-pass  filter,  and  outlining  a  step-by-step  design  procedure  drawn 
from  material  in  the  standard  text-books  on  the  subject  [3,4,5  J  . 

The  following  sections  describe  the  basic  theory  and  the 
design  procedure;  give  a  worked  example  and  discuss  some  practical 
aspects  of  the  working  filters. 

2.  BASIC  THEORY 

2 . 1  General 


Die  following  sections  are  a  summary  of  the  principles 
treated  in  greater  detail  in  references  4  and  5. 

In  this  discussion,  it  is  assisted  that  original  signal  of 
interest  is  an  analogue  signal  (e.g.  the  output  of  an  accelerometer) 
which  is  sampled  at  regular  intervals.  These  analogue  samples  are 
then  digitized  by  an  Analogue- to- Digital  Converter  (ADC) .  The 
data  to  be  digitally  filtered  are  the  successive  digital  words  from 
the  ADC  output. 

For  the  majority  of  applications,  the  unwanted  components 
of  an  analogue  signal  are  higher  in  frequency  than  the  desired 
components,  and  thus  a  low  pass  filter  is  the  required  format. 

Again,  a  flat  amplitude  response  in  the  pass-band  is  sore  important 
in  most  applications  than  a  linear  phase  response,  and  so  the 
Butterworth  type  of  filter  is  chosen. 

From  the  considerations  given  above,  the  digital  Butterworth 
low  pass  filter  will  fulfil  most  of  the  day-to-day  filtering 
requirements  and  thus  has  been  selected  as  the  design  example. 
However,  the  extension  of  the  procedures  to  other  types  of  filter 
(e.g.  chebychev  or  elliptical)  is  mot  difficult. 


A  digital  filter  is  a  digital  network  which  has  the 
property  of  being  frequency  selective.  It  works  by  multiplying 
current  (in  time)  and  delayed  samples  of  the  input  data  (and  in 
some  cases,  of  the  output  data)  by  coefficients  and  adding  the 
results  in  an  ordered  manner.  Filters  in  which  the  delayed 
samples  of  output  are  used  are  described  as  recursive  filters  and 
those  in  which  delayed  samples  of  the  input  are  used  are  known 
as  non-recursive. 

An  example  of  a  simple  recursive  digital  filter  is  shown 
in  FIG.  1.  The  input  sample,  X(N),  is  multiplied  a  coefficient 
bQ  and  added  to  the  preceding  input  san$>le,  X(N-l),  which  is 
itself  multiplied  by  a  coefficient  b1  to  give  the  intermediate 
sum: 

bQ  X (M)  +  bx  X(N-l) 

The  output  is  composed  of  the  intermediate  sum  added  to  the 
previous  output  sample  Y(N-l),  multiplied  by  the  coefficient,  a^. 
The  output  sample  Y(N),  is  thus  given  by: 

Y(N)  *  ax  Y(N-l)  +  bQ  X(iJ)  +  ^  X(N-l)  (1) 

2.2  The  Z-Transform 


A  mathematical  transform  which  has  been  found  to  be  useful 
in  the  representation  of  digital  networks  is  the  Z-transform.  The 
Z- transform,  X(Z)  of  a  sequence  of  N  values  of  X  (e.g,  the  samples 
of  a  transducer  output)  is  defined  by: 


X(Z)  *  l  X (N) Z 


where  Z  is  a  complex  variable.  If  Z  is  written  in  polar  form 
(Z  -  re^w,  where  w  is  angular  frequency)  and  if  r»l,  then  the  Z- 
transform: 

OD 

X(Z)  -  E  X(N)  e"1"1’ 

N--» 


is  the  sum  as  ths  discrete  Fourier  transfozsi  of  the  sequence. 

The  principal  merit  of  using  ths  Z-transform  in  the  synthesis 
of  digital  filters  is  that  the  Z-transform  of  s  delayed  sample 
X(N-l)  is  Z~1  timss  ths  z-transform  of  X(N).  Thus  ths  Z-transform 
snablss  digital  delays  to  be  represented  mathasMtically  in  a 
convenient  way.  For  example ,  FIG.  1  could  be  redrawn  with  the 
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unit  sample  boxes  replaced  by  Z*1  boxes  (FIG.  2).  If  FIG.  2  were 
much  more  conqplex  with  M  orders  of  delay  (i.e.  the  filter  order 
is  M)  instead  of  just  one  (e.g.  FIG.  3),  the  general  expression 
for  such  a  filter  would  be  4,5 

M  M 

Y(tf)  *  Z  a  Y(N-K)  +  Z  b  X(N-K)  (2) 

K*1  K  K«0  K 


If  M*l,  we  have  the  situation  in  FIG.  2. 


In  the  Z-domain,  the  input  and  output  of  a  digital  filter 
are  related  by  the  system  function,  H(Z)  which  is  defined  as: 


H(Z) 


Y  (Z) 
X(Z) 


where  Y(z)  is  the  Z-transform  of  Y (N)  and  X(Z)  is  the  Z- transform 
of  X(N). 


The  filter  will  in  general  satisfy  the  linear,  constant- 
coefficient  difference  equation 


M 


M 


Z  a  Y(N-K)  *  Z  b  X(N-K) 
K-0  K  K-0  K 


Taking  the  Z- transform  of  both  sides,  and  using  the  properties 
of  the  Z- transform,  we  obtain: 


M 

Z  a  Y  (N 
K-0  K 


,-K>] 


£  b  X(N 
K-0  K 


-K) 


m  r 


M 

£  a_  Z  Y(N-K)  «  £  b  Z  X(N-K) 

K-0  L.  — I  K-0  L  — 


M  M 

£  «_  Z"K  Y (Z)  -  £  br  Z"K  X(Z) 

K-0  K-0 


H(Z) 


Y(Z)  £  b  Z" 
K-0  K 


X(Z)  i  a  Z-K 
K-0  K 
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Thus  the  definition  of  a  digital  filter  is  a  matter  of 
determining  the  numbers  and  value  of  the  coefficients  and  a^. 
For  example ,  the  system  function  of  a  typical  sixth  order  digital 
filter  with  a  sanpling  frequency  of  60  Hz,  a  cut-off  frequency 
of  13  Hz,  and  0  dB  gain  the  pass  band,  is: 

H(Z)  *  0> 0147 (1+6Z~^-  +  15Z~2  +  20Z~3  +  15Z~4  +  6Z~5  +  Z~6 ) 
(1-0, +  1,0263Z~2  -  0,4243Z~3  +  0,182lz“4 

-  0.0305Z-5  +  0,0032Z-6) 

and  the  filter  is  completely  defined  by  this  expression. 

2.3  Flow  Graphs 

The  network  of  FIG.  4  can  be  more  simply  represented  by 
a  flow  graph,  as  in  FIG.  5.  Because  the  filter  network  is 
linear,  and  is  also  shift- invariant  [V]  ,  the  order  of  the  two 
ladder  networks  in  FIG.  4  can  be  reversed  without  changing  the 
input/output  relationship.  Reversing  the  ladder  networks ,  ^d 
eliminating  redundant  delays,  the  resulting  network  £fig.  5J 
has  the  minimum  number  of  delays  to  realise  the  general  filter 
and  is  called  a  CANONIC  form  of  the  filter. 

2.4  Filter  Design  Based  on  Analogue  Design  Methods 

Many  years  experience  in  the  design  of  analogue  filters 
has  resulted  in  well  tested  and  relatively  simple  design 
procedures  being  available*  It  has  been  found  convenient  there¬ 
fore  to  design  digital  filters  using  methods  based  on  the 
established  analogue  practice. 

In  analogue  filter  design,  the  system  function  which 
relates  the  output  to  the  input  is  normally  derived  in  the  S 
(or  Laplace)  domain,  and  is  of  the  general  form: 

M 

Y (S)  Ed*  SK 

K*0  * 

H(S)  -  -  -  -  (4) 

X(S)  E  C*  SK 

K«0 

where  Y (S)  is  the  Laplace  transform  of  the  output,  Y(t)  and 
where  X(S)  is  the  Laplace  transform  of  the  input,  X(t),  with  t 
being  the  time  variable,  and  M  being  the  filter  order. 


I 
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An  example  of  an  S-plane  system  function  is  the  fourth 
order  low-pass  filter  function  with  a  cut-off  frequency  of  6  Hz, 
and  unity  gain  in  the  pass* band  given  by: 

H(s)  =  _  0,1994 _ 

s4  +  1, 744S3  +  1, 5244S2  +  0,7819S +  0,2003 

The  required  digital  filter  system  function  is  the  equivalent 
analogue  filter  function  transformed  from  the  S-plane  into  the 
Z-plane . 


Several  methods  have  been  evolved  for  the  S-plane  to 
Z-plane  tranf ormation .  Each  of  the  methods  has  advantages,  but 
the  Bilinear  Transformation  method  has  been  selected  for  the 
following  examples  because  it  is  simple  and  straight  forward. 

2.5  Bilinear  Transformation 

Equation  3  [section  2.2]  is  a  differential  equation, 
and  the  bilinear  transformation  involves  integration  of  the 
differential  equation,  approximation  of  the  resulting  integral 
by  the  trapezoidal  rule,  Z- transforming  the  result,  and  solving 
for  H(Z)  [sj  .  By  inspection  of  the  expressions  for  H(S)  and 
H(Z),  which  result  from  the  calculation,  it  can  be  seen  that  the 
substitution  to  give  the  required  transformation  is: 

S  .  2-u-r.*>  (5) 

T  (l+Z”1) 

where  T  is  the  sample  period  of  the  data.  Substitutions  of  this 
form  are  generally  described  in  the  literature  as  bilinear 
transformations • 

In  the  calculations  which  follow,  the  data  points  are 
assumed  for  convenience  to  be  a  unit  time  apart.  Thus  the 
transformation  may  be  written: 


s  -  l  t1"2'1) 

(l+z-1) 


(6) 


The  effect  of  this  transformation  is  to  map  the  left  hand  side  of 
the  S-plane  into  a  unit  radius  circle  in  the  Z-plane  (FIG.  6), 
with  the  S-plane  imaginary  axis  being  transformed  into  the 
circumference  of  the  unit  circle.  The  frequency  response  of  the 
analogue  filter  is  evaluated  along  the  S-plane  imaginary  ( jfl) 
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axis,  while  the  frequency  response  of  the  digital  filter  is 
evaluated  along  the  unit  circle,  with  a  maximum  range  of  angle 
of  tir  radians.  In  practice,  this  range  of  ±tt  radians  corresponds 
to  ws^2f  where  ws  is  the  radian  sampling  frequency. 

A  problem  associated  with  this  type  of  transformation  is 
the  distorsion  in  frequency  as  we  go  from  the  S  to  the  Z  domain. 
Along  the  circumference  of  the  Z-plane  unit  circle,  Z  is  given 
by  e^w,  where  w  is  the  radian  frequency  variable.  From  equation  6, 


2(l-e~jw) 

(l+e-jw) 


Zj 


jfl 


FIG.  7  is  a  plot  of  the  digital  frequency,  v,  against  the  analogue 
frequency,  SI.  It  can  be  seen  that  the  closer  to  *rr  tnat  w  becomes, 
the  greater  the  difference  between  the  analogue  and  digital 
frequencies.  To  overcome  this  problem,  any  desired  digital  filter 
frequency,  w,  must  be  prewarped,  before  the  application  of  the 
bilinear  transformation,  to  obey  the  relation: 


SI  =  Z  tan 


(7) 


where  SI  is  the  prewarped  frequency  in  radians/sec. 

2.6  Analogue  Butterworth  Filters 

The  general  form  of  the  analogue  Butterworth  filter 
system  function  can  be  written  : 


H  (S) 


_ (nc)M _ 

(S-Pi)  (S-p2)  (S-P3)  . . .  (S-pj,) 


(8) 


where  Slc  is  the  cut-off  frequency  (in  radians/second),  defined 
as  the  frequency  at  which  the  system  function  is  3  decibels  (dP) 
less  than  the  pass-band  value  and  where  M  is  the  filter  order. 

The  values  of  p^  to  pM  are  the  roots  of  the  denominator  and  are 

the  poles  of  the  system  function.  When  plotted  in  the  S-plane, 

the  poles  lie  on  a  circle  of  radius  ftc,  spaced  ir/M  apart  on  the 

circumference,  ar\d  symmetrically  distributed  about  the  imaginary 
(jfl)  axis  [FIG.  8j  .  Thus  to  graphically  determine  the  right  hand 
side  of  equation  8,  all  that  is  required  is  a  knowledge  of  the 
filter  order ,  M,  and  the  cut-off  frequency  Slc . 
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In  general 
Butterworth  filter 


form,  the  frequency  response  of  an  analogue 
is  given  by  [sj 


|H(jQ)| 


2M 


where  Q  is  any  analogue  radian  frequency,  Qc  is  the  analogue 
radian  cut-off  frequency  and  M  is  the  filter  order.  The  filter 
altenuation  A  in  decibels  at  any  frequency  is  given  by  the 
relation  : 


A  =  20  log10  | H ( jft)  | 


=  20  log^Q  ^ 

/  i 

/ft 

V  nJ 

*  10  logio 

1 

(l+ft/ttc>2M 

Therefore,  10A/1°  =  1  + 

(ft/ftc)2M 

Recording  equation  9,  we 

obtain: 

ftc  =  ft  (lO7^10  -  l)"*5”  (10) 


If  our  filter  requirements  are  known  in  terms  of  the 
desired  altenuations  A^  and  A2  at  two  respective  frequency  points, 
and  we  can  determine  the  filter  order  from  two  simultaneous 
equations.  Eliminating  Qc,  we  obtain: 


M 


,  f  A2/10  Al/10 

login  I (10  ~  l)/(lo  -  1)J 
2  log^Q  (^2/^1^ 


(ID 


3.  DESIGN  STEPS 

The  tools  required  for  the  design  procedure  are  a  pocket 
calculator,  a  drawing  compass,  protractor,  linear  and  log-linear 
graph  paper.  The  method  involves  eight  steps  which  are: 

(i)  Specify  the  required  filter  by  selecting  two 
critical  analogue  frequency  points,  one  in 
the  pass  band  (w^)  and  one  in  the  stop-band  (w2) 


with  the  required  attenuation  at  each  of  these  points, 

(AX  &  A2) . 

if  the  sampling  frequency,  ws,  is  not  pre¬ 
determined,  one  should  be  chosen  such  that  it  is  higher 
than  the  stop-band  frequency,  . 

(ii)  Prewarp  the  critical  frequencies  (equation  7)  to  allow 
for  the  frequency  distorsion  inherent  in  the  bilinear 
transformation . 

(iii)  Determine  the  filter  order..  M,  from  eouation  11.  Round 
up  to  the  nearest  integer. 

(iv)  Determine  the  cut-off  frequency  from  equation  10. 

(v)  Graphically  determine  the  analogue  system  function, 

H (S) .  (See  section  2.6). 

(vi)  Apply  the  bilinear  transformation  (equation  6)  to 
H (S)  to  obtain  H(Z). 

(vii)  Draw  the  flow  graph  by  inspection  of  H(Z). 

Implement  the  flow  graph  as  a  computer  or  microprocessor 
program. 

It  should  be  noted  that  the  frequencies  w^  and  w2  in 
step  (i)  must  in  terms  of  the  digital  sampling  frequency,  w  ,  which  is 
assigned  the  value  of  2tt  radians/sec.  The  filter  characteristic  is 
periodic  in  2tt  ,  with  the  maximum  attenuation  being  at  tt  radians/sec 
(the  Nyquist  or  folding  frequency) .  The  filter  characteristic  above  tt 
will  be  a  mirror  image  of  the  characteristic  below  tt. 

4.  WORKED  EXAMPLE 


The  design  of  a  2-pole  digital  Butterworth  low-pass  filter 
is  used  to  illustrate  the  eight  step  procedure  outlined  above. 

STEP  1 


Assume  that  the  filter  requirements  are  a  pass  band 
attenuation  of  1  dB  at  1,5  Hertz  (Hz)  and  a  stop-band  attenuation  of 
15  dB  at  6  Hz.  Assume  a  sampling  frequency  of  60  Hz  giving  a  Nyquist 
frequency  of  30  Hz  which  is  equivalent  to  n  radians/ sec.  in  the 
digital  filter  frequency  domain.  The  critical  frequencies  must  be 
known  in  terms  of  the  sampling  frequency,  thus: 

wi  *  1*5  (2tt )  *  0,1571  radians/sec. 

60 


6,0  (2tt)  «■ 
60 


w  * 


0,6283  radians/ sec. 
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Thus  the  attenuations  =  1  dB  at  w^  -  0,1571  radians/sec. 

and  A2  =  5  dB  at  w2  '  0,6263  radians/sec. 

which  is  the  digital  filter  specification. 

STEP  2 

Prewarp  the  critical  frequencies 


'i 


=  2  tarn  0,1574  radians/sec. 
ft  2  "  2  tan^j?^  =  0,6498  radians/ sec. 
check  that  ft^  >  wl  and  ^2  >  w2* 


STEP  3 

Determine  the  filter  order. 

M  =  log10  [(100‘1-l)/g01,5-l) 

2  log10  (0,6498/0,15774) 

=  1,6831 

Round  M  to  next  hightest  integer  M  =  2. 

STEP  4 

Determine  the  cut-off  frequency 

15  -V 
ft  *  ft-  (10  '  -  1) 

c  2 

*  0,2762  radians/sec. 

check  that  this  is  close  to  the  requirements. 


STEP  5 


Determine  the  S-plane  (Laplace)  transfer  function  from 
S-plane  pole  plot  for  a  two-pole  analogue  Butterworth  filter 
FIG.  9j  . 


To  determine  the  pole  values,  draw  a  circle  of  radius 
0,2762  units  and  mark  the  poles  at  it/  apart  on  the  left-hand  side 
of  the  circle  only.  The  poles  must  be  symmetrically  placed  about 


i 
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the  imaginary  (vertical)  axis  [see  FIG.  8^]  .  The  pole  positions  are 

then  measured  from  the  diagram  and  the  S-plane  system  function  is 
obtained. 


If  greater  accuracy  is  required,  the  pole  positions  may 
be  determined  by  calculating  the  projection  of  the  pole-to-origin 
distance  on  the  real  and  imaginary  axes  since  the  radius  (0,2762)  and 
the  angles  are  known  precisely. 

The  system  function  (by  graphical  methods)  is: 


H  (S) 


_ (0, 2762) 2 _ 

(S+0, 195+j0, 195)  (S+0, 195- jO, 195) 


=  0>0763 

S  +0,395+0,0761 

STEP  6 

Apply  bilinear  transformation. 

Substitute  2(1-Z~1)  for  S  in  the  denominator 
1+Z-1 

of  the  expression  for  H(S)  above 
S2  +  0,395  +  0,0761 

=  4 (1-2Z~1+Z-2)  +0, 78-0, 78Z~2+0, 0761+0, 1522Z_1+0, 0761Z"2 

r  -  .....  ...  -  — - 

1+2Z~1+Z~2 

-  4, 8561(1-1, 6161Z~1+0, 6788  Z~2j 
1+2Z-1+Z-2 

Therefore  the  Z-plane  system  function  becomes 

0,0157  (1+2Z"1+Z-2) 

H  (Z)  *=  - * - 5 - ~ 

1+1 , 6161Z” 1-0, 6788Z- 2 

STEP  7 

From  the  system  function,  draw  the  flow  diagram  (fig.  1(3  . 
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STEP  8 

Implement  the  filter  as  a  computer  program.  FIG.  11  shows 
the  results  of  implementing  the  filter  on  a  large  digital  computer  and 
with  a  microprocessor.  As  can  be  seen  from  the  plot,  the  objectives, 
as  defined  in  Step  1  have  been  achieved. 

The  further  worked  examples  of  digital  Butterworth  filters 
fife  ffivered^n  gffnma£v  in  Apjjendices^l  and  2^. 

5,  DISCUSSION  OF  RESULTS 


In  implementing  the  Z~plane  system  function,  H(Z),  care 
must  be  taken  to  ensure  the  accuracy  of  the  denominator  coefficients. 

For  example  [FIG.  12]  ,  a  5%  error  in  one  of  the  denominator  coefficients 
of  a  4-pole  filter  produces  a  complete  breakdown  of  the  filter.  A 
1%  error  in  one  denominator  coefficient  produces  a  significant  error. 

It  is  therefore  desirable  to  aim  at  an  accuracy  of  £0.1%  in  the 
denominator  coefficients.  The  denominator  coefficient  sensitivity 
can  be  reduced  somewhat  by  choosing  a  less  efficient  structure  13,4] 
to  the  canonic  form  used  in  the  example  above. 

Related  to  the  problem  of  coefficient  accuracy  is  the 
question  of  rounding  errors  due  to  finite  register  length  in  micro¬ 
processor  implementations  of  digital  filters.  For  many  applications, 
the  calculation  speed  of  the  most  commonly  available  microprocessors 
is  not  sufficient  to  allow  software  floating  point  arithmetic  to  be 
used,  and  so  integer  arithmetic  must  be  used  instead.  The  errors 
associated  with  the  finite  (8  or  16-bit)  register  length  must  be 
evaluated  for  each  particular  filter  and  the  general  methods  are 
discussed  in  references  3  and  4. 

An  alternative  approach  is  to  use  one  of  the  special 
arithmetic  chips  (e.g.  advanced  Micro  Devices  AM9511  or  similar)  to 
perform  the  floating  point  operations  or  the  integer  multiply  and 
divide  functions. 

A  minor  point  which  has  been  found  useful  in  the  micro¬ 
processor  applications  is  to  scale  the  coefficients  to  integer  values 
by  multiplying  by  numbers  which  are  some  power  of  2  rather  than  by 
some  power  of  10.  For  example,  the  coefficient  1,616  could  be 
converted  to  an  integer  value  by  multiplying  by  1000.  However,  a 
result  of  sufficient  accuracy  can  be  obtained  by  multiplying  by 
1024  (*2l0)  which  involves  a  10-bit  shift  left  rather  than  a  complex 
multiplication  by  1000. 
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6.  CONCLUSION 

The  majority  of  practical  filtering  tasks  associated  with 
data  gathering  can  be  accomplished  using  low-pass  filters,  and  this 
memorandum  provides  a  simplified  suitmary  of  one  of  the  several  text¬ 
book  approaches  to  digital  filter  design,  forked  examples  are  included 
to  illustrate  and  amplify  the  summary.  The  extension  of  the  low-pass 
methods  to  cover  the  high  pass  and  band-pass  cases,  as  well  as  filter 
types  other  than  Butterworth,  is  discussed  in  a  straight  forward 
manner  in  the  referenced  texts. 
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APPENDIX  1 


SUWiARY  OF  THE  DESIGN  OF  A  4-POLr 
DIGITAL  BUTTE RWORTH  LOW-PASS  FILTER 


The  sampling  frequency  is  assumed  to  be  60  Hz,  equivalent  to 
2n  radians/sec.  in  the  digital  frequency  domain.  The  attenuation 
requirements  are: 


0,5  dB  at  4  Hz;  -27  dB  at  12  Hz. 


STEP  1 


A1  «  0,5  dB  at  wj  =  0,4199  radians/ sec. 
A2  m  dB  at  Wj  *  1,2566  radians/sec. 

STEP  2 

*  0,4251  radians/sec. 

«  1,4530  radians/sec. 


STEP  3 


M  -  3,384 

Round  to  next  heichtest  integer  m  *  4. 


STEP  4 


flc  *  0,6632  radians/sec. 


STEP  5 


STEP  6 


STEP  7 


H  (S) 


_ 0,1994 _ 

S4+l, 744S3+1, 5244S2+0, 7819S+0, 2003 


H(Z)  - 


(0,005273) (1+4Z~1+6Z~2+4Z~3+Z“4> _ 

1- 2 , 3265Z  *+2 , 2480Z" 2- 1 , 0161Z" 3+0 , 1794Z*4 


Plow  graph  appears  in  FIG.  13. 

STEP  8 


Results  of  computer  implementation  appear  in  fig.  14. 
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APPENDIX  2 

SUMMARY  OF  THE  DESIGN  OF  A 
6-POLE  DIGITAL  BUTTERS  10 RTH  LOS'?- PASS  FILTER 


The  sampling  frequency  is  assumed  to  be  60  Hz  (=  2*  radians/ 
sec.  in  the  digital  frequency  domain) .  The  attenuation  requirements 
cure  : 


0,2  dB  at  9  Hz  and  40  dB  at  20  Hz. 

STEP  1 

A^  *  0,2  dB  at  *  0,9425  radians/sec. 
A2  *  40  dB  at  w2  *  2,0944  radians/sec. 

STEP  2 

*  1,0191  radians/sec. 

=*  3,4641  radians/ sec. 


STEP  3 

M  -  5,0122,  rounded  up  to  6. 


I  STEP  4 

*  1,6079  radians/ sec. 


STEP  5 

17,2804 

H  (S)  - 

S6+6, 21S5+19 , 2965S4+38 , 0149S3 


+49, 9034S2+41,4914S+17, 2387 

STEP  6 

H(2>  m  0,014701 (l+6Z"1+lSZ~2+202~3+lSZ~4+6Z~5+Z~6)  - 

.1-0,8182Z_1+1,0263Z~2-0,4243Z~3+ 


APPENDIX  2  (CONTD.) 


STEP  7 


Flow  graph 


STEP  8 


Implement  on  computer.  Frequency  characteristic  of  the 
resultant  filter  is  shown  in  FIG.  16. 


DIRECT  FORM  REALIZATION  OF  GENERAL  DIGITAL  FILTER  OF 
ORDER  M. 


GRAPHACAL  REPRESENTATION  OF  BILINEAR  TRAILS  FORM  AT  I  ON' 


FREQUENCY  DISTORTION  DUE  TO  BIT. T INEAR  TRANSFORMATION 


I  MAG 

(jB) 


Pole  positions 


PIG.  8  'BUTTERWORTH  CIRCLE',  TYPICAL 

EXAMPLE.  POLE  PLOT  OP  4 -POLE  BUTTERWORTH 

ANALOGUE  LOW  PASS  FILTER,  n  -  0,5  RADS/SEC. 

c 


a1  =  1,6161  bQ  =  0,0157 

a  =  -0,6788  b,  =  2,00 

b2  »  1,00 


PIG.  10  PLOW  DIAGRAM  FOR  2-POLE  DIGITAL  BUTTERWORTH  LOW  PASS 
FILTER 
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FIG.  13  FLOW  GRAPH  OF  A  4-POLE  DIGITAL  BUTTERWORTH  LOW-PASS  FILTER 

(DIRECT  FORM  II  (REF.  5)) 


Gain  (dB) 


FIG.  14  FREQUENCY  CHARACTERISTIC  OF  4-POLE  DIGITAL 
BUTTERWORTH  FILTER. 
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